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Abstract. We present a study of the hard X-ray spectrum ( > 15 keV) of different classes of Seyfert galaxies observed 
with BeppoSAX/PDS. Using hard X-ray data, we avoid absorption effects modifying the Seyfert emission and 
have direct access to the central engine of these sources. The aim of this study is first to characterize the general 
properties of the hard X-ray spectrum of Seyfert 1, 1.5 and 2 galaxies and secondly to compare their intrinsic 
emission to test unified models according to which all the classes have the same nucleus. 

We compute the average spectrum of 14 Sy 1, 9 Sy 1.5 and 22 Sy 2 galaxies observed by the PDS (15-136 keV). 
The average spectrum of Sy 1 differs from that of Sy 2, the first requiring the presence of a high energy cutoff 
which is absent in the second. We also show that the reflection component is possibly more important in the Sy 
2 emission. The nature of Sy 1.5 galaxies is ambiguous as they have a negligible reflection component (like Sy 1) 
and do not require a cutoff (like Sy 2). 
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1. Introduction 

According to unified models l|Antonucci fc Milled. Il985(l . 
Seyfert 2 galaxies harbor a bright Seyfert 1 nucleus which 
| is hidden from our view by an optically and geometri- 
cally thick obscuring torus. Consequently, the difference 
between type 1 and 2 Seyfert galaxies would be only due 
to viewing angle. Several observational facts in favor of 
unified models have been found: 



— detection of Polarized Bro ad Lines (PBL hereafter) in 



several Seyfert 2 galaxies (Antonucci & Milled. Il985 


Miller & Goodrich! 


199(1 iTrar, 


119951 lYountr et al. 


199rl lAwaki et a 


, EpOO; 


Moran et all I200C 


Alexandel l200ll 


Gu et all l200ll iLumsden et al. 


2001), interpreted as scattering 


; of the Broad Lines 



Region emission by warm material placed above an 
absorption torus. 



— large absorbing column densities obse rved in the hard 
X-ray spectra of Seyfert 2 galaxies ijGuainazzi et alJ 
(200l| and references therein). 

— suspicion of anisotropic ionization radiation from the 
nucleus through the structures in the light of [OIII] line 

Send offprint requests to: S. Deluit, 
e-mail: Sandr ine . DeluitOobs . unige . ch 



at 5007 A (e.g. in NGC 5252 l|Tadhunter fc Tsvetanovl 

EH). 

X-ray data have played a major role in support of uni- 
fied models. The X-ray spectra of Seyfert galaxies are 
known relatively well and are composed at energies above 
3 keV of a pow er law, a Fe line and a Compton reflec- 
tion c o mponent ijNandra fc Pounds! lll994h.lGondek et alJ 
lll996l) iNandra etalJ l|l997|h IZdziarski et all (|l999|h 
Zdziars ki et al.l l|2000|) ). Some questions remain like the 
shape and the presence of a cutoff at higher energies. 
Furthermore, some observations are inconsistent with uni- 
fied models at least in their simplest formulation. High 
energy spectra of a sample of Seyfert 2 galaxies observed 
by OSSE suggest an intrinsic difference in the hard X- 
ray emission of both classes of Seyfert galaxies. The two 
classes differ in the slope of the intrinsic power law, in the 
energy of the exponential cutoff, in the amount of absorp- 
tion and reflection (Zdziarski, 1995). 
Hard X-ray spectra are a powerful tool to probe and test 
unification models by providing a direct view on the cen- 
tral engine emission mostly free of absorption effects. The 
form of the individual hard X-ray spectra is poorly con- 
strained due to limited photon statistics. This can be im- 
proved by considering the average spectrum of a sample 
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to extract the properties of a class. 

We study the hard X-ray emission (above 15 keV) of 
Seyfert 1, 1.5 and 2 galaxies in order to investigate several 
questions: 

— Are the high energy index slope of Seyfert 2 and 
Seyfert 1 identical? 

— Do the Seyfert 1 and 2 spectra have a cutoff? In ther- 
mal Compton models the high energy cutoff provides 
a reliable estimate of the comptonizin g material tem- 
perature l|Rvbicki k, Lightmanl 1197^1 . Therefore, dif- 
ferences in cutoffs or their absence indicate different 
comptonizing medium temperatures. 

— What is the importance of Compton reprocessed emis- 
sion for each class? Could the reflection be the only 
cause to the observed difference between Seyfert 1 and 
Seyfert 2? 

To answer these questions, we compute the average PDS 
spectrum of each class. We consider for inclusion in the 
Seyfert 2 sample only galaxies for which the absorption 
column density is less than 7T0 23 cm~ 2 to ensure that, 
above 15 keV, absorption is less than 50% of the flux. 
We are thus in the position to characterize the intrinsic 
emission of the various types of objects. 



such objects from consideration. In summary, we consider 
only objects for which c^ c ' tNh x e - rJ ph{E)N H j s l ar g e r 
than 0.5 at 15 keV, ctt and a p h are respectively the 
Thomson and the photoelectric cross section. This means 
that we keep only objects having a column density less 
than 7-10 23 cm- 2 . 

The role of absorption in our Seyfert 2 class study is ex- 
haustively explained in the appendix B. 
We present in Figure 1 the column density distribution for 
the Seyfert 2 galaxies sample. The majority of the objects 
have Njj between 10 22 cm -2 and 10 23 cm~ 2 . 

Distribution of Seyfert 2 hydrogen column densities 



2. Selection procedure of our sample 



We start by including all (i.e. 68) Seyfert PDS spectra 
available in the public archive of BeppoSAX Data Center 1 . 
We have investigated the possibility of source confusion 
in the PDS aperture by looking in the MECS images and 
excluded two confused objects from this study. We then 
form three subsamples characterizing Sy 1, Sy 1.5 and Sy 
2 galaxies respectively. 



Log <N H ) 

Fig. 1. Distribution of the hydrogen column density of 
sources composing our Seyfert 2 sample 

We present in Figure 2 the hydrogen column density 
distribution versus the redshift for Seyfert 2 galaxies. 



2.1. The nature of objects composing the subsamples 

We follow the optical classification of Seyfert galaxies indi- 
cated in the NASA/IPAC Extragalactic Database 2 (NED) 
and exclude LINERs for the ambiguity of their classifica- 
tion. We include Sy 1 and Sy 1.2 galaxies in the Seyfert 1 
subsample, Sy 1.5 in the Seyfert 1.5 subsample and Sy 1.9 
and Sy 2 in the Seyfert 2 subsample. We exclude Seyfert 
1.8 for which the classification is too ambiguous between 
Seyfert 1.5 and Seyfert 2 galaxies. 



2.2. Absorption effects in Seyfert 2 galaxies emission 

Hard X-rays give a direct access to the intrinsic emis- 
sion provided that we avoid absorption effects which con- 
taminate the Seyfert 2 galaxies spectra at lower ener- 
gies. Objects with Njj>1.5T0 24 cm -2 (so called Compton 
thick) have a Thomson optical depth larger than one. 
Electron scattering modifies their emission even in the 
hard X-ray energy domain. We therefore eliminated all 



Distribution of the Seyfert 2 column density along the Redshift 



1 http:/ /www. asdc.asi.it/bepposax/ 

2 http:/ /nedwww. ipac.caltech.edu/ 
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Fig. 2. Distribution of Seyfert 2 hydrogen column density 
along the redshift 

We note a large number of high column densities at low 
redshift. This is probably an observational artefact due 
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to the fact that distant and highly absorbed Seyfert 2 
galaxies are hardly detectable. 

2.3. Data quality criterion 

The third criterion we introduce is the quality of the data. 
We calculate the integrated signal to noise ratio of each 
object in the PDS energy range (15-136 keV). The data 
are used only when the integrated signal to noise ratio is 
equal to or larger than 2o\ This threshold is high enough to 
conclude to a detection by the PDS detector for a known 
source. 

3. Properties of our samples 

3.1. Redshift distributions 

We present in Figure 3 the redshift distribution for the 
three classes of Seyfert galaxies. 



Distribution of the Sy 1 integrated Signal to Noise ratio in the PDS energy range (1 5—1 36 keV) 
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Fig. 3. Redshift distribution for the three classes of 
Seyfert galaxies 

Most of the objects have a redshift less than 0.05 (93% 
for Sy 1, 89% for Sy 1.5 and 82% for Sy 2 galaxies). The 
difference between the redshift distributions of Sy 1 and 
Sy 2 galaxies may be due to biases in object selections 
when planning X-ray observation programmes. Indeed, 
absorbed objects being more difficult to observe will be 
looked for at small distances. 



3.2. Comparison of the data quality 

We determine the integrated signal to noise ratio in the 
PDS energy range (15-136 keV) for each object (Figure 
4). 

The majority of the objects have an integrated signal 
to noise ratio less than 30cr (79% for Sy 1, 78% for Sy 



Fig. 4. Distribution of the integrated signal to noise ratio 
for the three classes of Seyfert in the PDS energy range 



1.5 and 82% for Sy 2 galaxies). The Seyfert 1, 1.5 and 2 
samples are dominated by low signal to noise spectra. 

3.3. Flux distributions 

We measure the integrated flux in the PDS energy range 
of each object (Figure 5). The majority of the sources have 
a hard X-ray flux less than 15-10 _11 ergs _1 cm~ 2 between 
15 keV and 136 keV. 



Distribution of the Sy 1 Flux in the PDS energy range (15-136 keV) 
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Fig. 5. Flux distribution for the three classes of Seyfert in 
the PDS energy range 
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The proportion of objects with a flux less than or equal 
to 15-10- 11 ergs- 1 cm- 2 is of 93% for Sy 1, 100% for Sy 
1.5 and of 82% for Sy 2 galaxies. 

3.4. Luminosity distributions 

From the flux we calculate the integrated luminosity using 
Ho=75 kms -1 Mpc -1 and qo=0.5 as cosmological param- 
eters. The results are shown in Figure 6. 



Distribution of the Sy 1 Luminosity in the PDS energy range (15-136 keV) 
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Fig. 6. Luminosity distribution for the three classes 
Seyfert in the PDS energy range 
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The majority of the objects have a hard X-ray lumi- 
nosity less than or equal to 30T0 43 ergs _1 between 15 keV 
and 136 keV (93% for Sy 1, 100% for Sy 1.5 and 91% 
for Sy 2). In the Seyfert 2 subsample, two objects have a 
luminosity above SO-lO^ergs" 1 : IRAS F12072-0444 with 
a luminosity of ei^-lO^ergs" 1 and IRAS 23060+0505, 
with a luminosity of lll-10 43 ergs _1 . These two galax- 
ies are warm Ultra Luminous Infrared Galaxies and are 
probably quasars regarding their hard X-ray luminosity. 
However, they are classified as Seyfert 2 galaxies in NED 
and therefore kept in our sample. 

3.5. Excluded objects 

Several objects have been excluded from our initial sam- 
ple. We review below all the dismissed sources and the 
reason of their exclusion. 

1. NGC 4151: with an integrated signal to noise ratio of 
227.5cr and an integrated flux in the PDS energy range 
of 58-10 _11 ergs _1 cm -2 , this object would strongly 
dominate the sample. We will discuss this bright 
Seyfert 1.5 galaxy in the paper II of this study where 



indiv idual objects are discussed ijDeluit fc Courvoisieri 
12001 . 

2. Mrk 1073: this object is not detected in the MEGS 
and LECS instruments and has a very high signal in 
the PDS. We excluded it from our sample because its 
PDS spectrum is probably in part due to a nearby 
source. 

3. NGC 1386: a large fr action of the PDS flux is due to 
the nearby NGC 1365 l|Maiolino et allll998l) . 

4. MCG 5-18-002: we excluded this source in reason of 
its uncertain absorpt ion given as Njj >10 25 cm~ 2 by 
iMaiolino et al.l l)l998|) . The value is based only on the 
[OIII] flux and requires further confirmation. 

5. Considering the data quality criterion, we also exclude 
four Seyfert 1 (TON S180, GQCOMAE, IRAS 
13224-3809 and Mrk 478) and four Seyfert 2 (NGC 
3081, IRAS 11058-1131, IRAS 10126+7339 and 
NGC 4941). All these objects have an integrated sig- 
nal to noise ratio less than 2a. 

3.6. The final sample 

The initial sample of all Seyfert galaxies observed with 
BcppoSAX included 68 objects among 18 Sy 1, 10 Sy 1.5 
and 40 Sy 2 galaxies. 

Our sample is finally composed of 67 observations for a to- 
tal of 45 Seyfert galaxies of which 14 are Seyfert 1 galaxies, 
9 are Seyfert 1.5 and 22 are Seyfert 2 galaxies. 
The list of objects composing our subsamples, their char- 
acteristics and their spectral properties is given in Table 
A.l. for Sy 1 and Sy 1.5 galaxies samples and in Table 
A. 2. for the Sy 2 galaxies sample. 

4. Data Analysis 

The BeppoSAX instruments cover a wide spectral band 
and consist of four co-aligned Narrow Field Instruments 
(NFI) plus two Wide Field Cameras perpendicular to the 
axes of the NFI and looking in opposite directions. 
The NFI is composed b y a Low Ene rg y Con centrator 
Spectrometer (LECS: IParmar et all (|l997h L three 
Medium Energy Concentrator Spectrometer (MECS: 
iBoella et all Jl997t>1 and a Phoswich Detector Counter 
CPDS: iFrontera et alJ l|l997|) L LECS and MECS have 
imaging capabilities and operate respectively in the 0.1-5 
keV and 2.0-10 keV. The instrument used in this study, 
the PDS, is an instrument using rocking collimators, the 
energy band is between 15-200 keV. The PDS consists 
in four phoswich units operating in collimator mode 
with two of them pointing towards the source, the other 
two away. The two pairs switch on and off source and 
the net source spectra are obtained by subtracting the 
'off (which represent the background) from the "on" 
counts. The PDS data available on the web are already 
background subtracted. A more detailed description 
of the background and of the subtraction procedure is 
available from http://www.sdc.asi.it/software. 
The data analysis has been performed with the XSPEC 
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ijArnaudl Il996l) version 11 package and using the last 
PDS response matrices released by the SAX Data Center. 
We kept in the input and output files the same binning 
for the average spectrum as was available for the original 
data of each source. For a source observed several times, 
we first compute its average spectrum. We then compute 
the average subclass spectrum taking into account of the 
exposure time of each data set. The energy band is from 
15 to 136 keV. 

5. Spectral properties of the different classes of 
Seyfert galaxies 

We present in Figure 7 the average counts spectrum of the 
three classes of Seyfert galaxies. 



Average Spectra of the three classes of Seyfert Galaxies 




20 50 100 

channel energy (keV) 



Fig. 7. Average spectra of the three classes of Seyfert 
galaxies: Sy 1 (triangle), Sy 1.5 (square) and Sy 2 galaxies 
(circle) 

The spectrum of each class has an integrated signal to 
noise ratio of 96a for Seyfert 1 galaxies, 74.4cr for Seyfert 
1.5 galaxies and 94a for Seyfert 2 galaxies. No normaliza- 
tion has been applied to obtain the average spectra shown 
in this figure. 

We compare the Seyfert classes following two ways. First 
through a model independent method comparing directly 
the counts spectra. Secondly through a model dependent 
method using a fitting procedure to investigate the phys- 
ical processes at the origin of the emission. 

5.1. Model Independent Method: comparison of the 
normalized counts spectra 

We compare the different average spectra normalizing 
them at 15 keV. This allows us to compare their emission 
independently of the average flux. The result is shown in 
Figure 8. On the left column, we present the normalized 
counts spectra of the three pairs of classes, on the right 
side we give the ratio of each pair. 



5.1.1. Comparison between Seyfert 1 and Seyfert 2 
emission 

The two spectra are similar below 40 keV confirming that 
the Sy 2 galaxies emission is not significantly modified by 
absorption effects. Above 40 keV, the Sy 2 emission domi- 
nates that of Sy 1. The difference between the two classes 
increases above 90 keV. Indeed, the hard photons popu- 
lation is more important in Sy 2 than in Sy 1 galaxies 
showing that the observed Sy 2 spectrum is harder than 
Sy 1. The deficit of hard photons in the Sy 1 class can be 
explained either by the presence of a cutoff in its emission 
and/or by a reflection process occurring in the Sy 2 emis- 
sion. The possibility that the Sy 2 intrinsic power law is 
harder than Sy 1 cannot be excluded. The Sy 1 average 
emission is 42% weaker than that of Sy 2 around 100 keV 
where they differ most in the normalized ratio. 



5.1.2. Comparison between Seyfert 1 and Seyfert 1.5 
emission 

The spectra are identical below 40 keV. The Sy 1.5 emis- 
sion slightly dominates between 40 and 65 keV. At higher 
energies, the Sy 1.5 spectrum contains a more important 
hard photons population. As in Sy 2 galaxies, the Sy 1.5 
emission domination above 80 keV can be produced by a 
cutoff in Sy 1 or by a harder Sy 1.5 primary power law. 
The Sy 1 average emission is 41% weaker than that of Sy 
1.5 at ^100 keV where they differ most in the normalized 
ratio. 



5.1.3. Comparison between Seyfert 1.5 and Seyfert 2 
emission 

The two spectra remain very close along the PDS en- 
ergy band. Indeed, the maximal deviation between the 
two spectra in the normalized ratio does not exceed 20% 
at ~40 keV where the Sy 2 average emission is marginally 
higher than that of Sy 1.5. Thus, we consider in a first 
approach that the spectrum of Sy 1.5 and Sy 2 present 
notable similarities. 



5.1.4. Summary 

The main result of this model independent study is that 
Sy 2 differ from Sy 1 galaxies and that Sy 1.5 and Sy 2 
are similar. In particular in the high energy domain where 
the Sy 1.5 and Sy 2 classes present a more important hard 
photons population compared with Sy 1 galaxies. The pro- 
cesses that can be responsible for the observed differences 
between the three classes is investigated in the next sub- 
section. 



6 



Sandrine Deluit & Thierry Courvoisier: The intrinsic emission of Seyfert galaxies 




1 cf 4 L- ■ i ■ ■ i 1 0.4 L . . , . ._ 

20 40 60 80 100 120 20 40 60 80 100 120 



Energy (keV) Energy (keV) 

^ . . . 1 1.4 r- ' . ■ 




Energy (keV) Energy (keV) 



Fig. 8. Left: three pairs of Spectral Energy Distributions normalized at 15 keV. Right: ratios of these SEDs. 



Table 1. Best fit parameters for simple models 



Model 


Seyfert 


r 


Inbreak / cuto f f 


r 2 


xVd.o.f. 


xi 


FTEST 




Class 


Index Slope 


(keV) 


Index Slope 




reduced 


Probability value 


Power Law 


1 

1.5 
2 


2.03 

l qn+ 05 
r.»u_ 05 

1.90 






50.0/12 
16.7/12 
44.2/12 


4.17 
1.39 
3.68 




Broken Power Law 


1 


i.oo_q 2i 


OB+V.U 
28-3.8 


z - ZJ -0.09 


14.7/10 


1.47 


2.2-10 -3 




1.5 

2 


1 62+ 021 

± -"^-0 .38 


zo -5.4 
n 4 +4.3 


2 01+° 13 
2 04+ - 07 


10.3/10 
13.9/10 


1.03 
1.39 


5.0- 10- 2 

3.1- 10 -3 


Cutoff Power Law 


1 




70+f 7 




15.2/11 


1.38 


4.5-10 -4 




1.5 


1 64+ 011 

J -°*— 0.19 






13.3/11 


1.21 


5.8-10~ 2 




2 


1 4^+0-15 


87±g 




17.6/11 


1.60 


1.8-10 -3 



Note: 

The uncertainties correspond to 90% confidence level based on a Ax 2 =2.7 criterion llLampton et alj . ll976l) 



5.2. Model Dependent Method: investigating the 
physical processes into Seyfert galaxies engine 

To interpret our data in terms of models taking the PDS 
spectral response matrices into account, we fit the spectra 
with different theoretical models. We first fit the spectra 
with simple models (power law, broken power law and 
cutoff power law). We then consider reflection with the 
PEXRAV model described below. 
All the results are summarized in Tables 1 and 2. 

5.2.1. Fitting of the spectra with simple theoretical 
models 

The power law model gives the worst statistical results for 
all the classes. For Sy 1 and 2 galaxies, this model is sta- 



tistically unacceptable with a reduced x 2 greater than 2 
indicating that their spectra require an additional compo- 
nent. For the Sy 1.5 galaxies, this model is more adapted. 
Sy 1 emission requires a cutoff considering both the statis- 
tical result of the \ 2 test and the FTEST probability, con- 
firming our model independent analysis. The presence of 
a cutoff is unprobable in Sy 1.5 and Sy 2 galaxies. Indeed, 
the x 2 method gives the best result for a broken power law 
(x 2 =L03 for Sy 1.5 and x 2 =L30 for Sy 2) compared with 
the cutoff model (x 2 =L21 and x 2 = 1-60 respectively). 
In subsection 5.1, we noted two maxima in the ratio be- 
tween Sy 1 and Sy 2 spectra around 40 keV and 100 keV. 
These features may be due to a combination of a cutoff in 
Sy 1 (here found at 70j^fy keV) and/or a reflection hump 
in Sy 2 class that we investigate in the next section. 
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PEXRAV PARAMETERS 



Seyfert 


r 


Efolded 


R 


xVd.o.f. 


xi 


cos 9 


Class 


Index Slope 


(keV) 


Reflection 




reduced 


fixed 


1 

1.5 


1 92 +01J<i 

0.49 

1 9fi+ 006 


9.8-10*] 


0.93+^° 


13.2/10 
10.8/10 


1.32 
1.08 


0.45 
0.45 


2 


1 qq+o.05 

i,y3 -0 45 


2-4-10f ol 


1.4715;?? 


13.1/10 


1.31 


0.45 


1 

1.5 


-i oq+U.28 


nni+NC 
zzi -158 

8-2-10f a] 


0.48+-- 

f) QO+0.38 


13.1/10 
10.7/10 


1.31 
1.07 


0.87 
0.87 


2 


2 oo+ 005 


3.6-10f o] 


9 oq+1.96 
z.oo_ 49 


13.0/10 


1.30 


0.30 


Seyfert 


r 


Efolded 


R 


xVd.o.f. 


xi 


cos 9 


Class 


Index Slope 


(keV) 


Reflection 




reduced 


free 


1 

1.5 


1 - yl -0.47 
1 Q5+009 


l.(H0? o] 


53 +NC 

U.OO_ 27 

n 36 +NC 

U.OD_ 23 


13.2/9 
10.7/9 


1.47 
1.19 


0.81 
0.95 


2 


2 oo+ 005 


2.8-10f o] 


2.55„ 2 .i3 


13.0/9 


1.44 


0.28 



Notes: 

The uncertainties correspond to 90% confidence level based on a Ax 2 =2.7 criterion llLampton et alJ . ll976t) . 
NC indicates that no constraint has been found on this limit. 

[a] the value indicates that the high energy cutoff is far above values that can be measured by the PDS showing the absence of 
a measured cutoff. 



5.2.2. Fitting of the spectra with the reflection model 
PEXRAV 

Previously, we suspected the presence of a reflection 
component in the Sy 1.5 and Sy 2 spectrum. We in- 
vestigate this hypothesis us ing the PEXRAV model 
ifMaedziarz k Zdziarskilll995|) in XSPEC. This model cal- 
culates the expected X-ray spectrum when a point source 
of X-rays is incident on optically thick, mainly neutral 
material. The parameter R measures the reflection com- 
ponent. If the primary source is isotropic and neither the 
source nor the reflector are obscured, then R can be linked 
to the solid angle seen through the reflector as R~^. This 
model and in particular the reflection component arc in- 
clination dependent. We apply this model in two ways: 
first, we fix the inclination from the normal of the disk at 
a medium angle to test the occurrence of the reflection as 
well as the presence of other components independently 
to any strong inclination influence. We then change the 
angle value fixing or letting it free in order to understand 
its role and to obtain the inclination value for which we 
have the best fit. The results are shown in Table 2. 

1. Intermediate inclination analysis 

The inclination is initially fixed at an intermediate 
value of cos #=0.45. 

The x 2 test indicates that PEXRAV significantly im- 
proves the data fitting of Sy 1 and Sy 2 galaxies com- 
pared with simple models, not however for the Sy 1.5 
spectrum that remains best represented by a broken 
power law. Therefore, reflection may play a role for 
the two main Seyfert classes and be absent in Sy 1.5 
emission. This is confirmed by the reflection amount 
which can be non existent for Sy 1.5 (R=0.73_o.7i) 
whereas its presence is strongly suggested in Sy 2 emis- 
sion (R=1.47_o.77) at the 2a level. The primary power 



law indices of the three classes are within 1<7, however 
Sy 2 emission could be softer compared with the Sy 1 
index slope ^=1.991.^5 for S Y 2 against r=1.92±g^j 
for Sy 1). 

A strong difference appears in the high energy cutoff 
which is relatively well identified for Sy 1 (E c =238_i76 
keV) and totally absent for Sy 1.5 and Sy 2 classes. 
The main result is that the spectrum of Sy 2 galaxies 
differs from that of Sy 1 with the absence of a cut- 
off and a possible presence of an important reflection 
component in the Sy 2 emission. 
2. Inclination influence 

We change the inclination value to test unified models 
which suggest that Sy 2 galaxies are seen edge-on. 

(a) Fixed inclination at different angles 

We fix the inclination at cos 0=0.30 and cos 0=0.87. 
The influence of the inclination parameter is not 
crucial on the final result of the \ 2 test. The best 
fit for Sy 1 and Sy 1.5 is obtained for cos 0=0.87, 
an angle of 30° from the normal of the disk. The 
reflection component value is lower for cos 0=0.87 
than for an inclination of cos 0=0.45 as predicted 
by the fact that the reflection is inclination depen- 
dent in the PEXRAV model. We also confirm that 
we obtain for the Sy 2 galaxies the best % 2 value 
for cos 0=0.30 namely an angle of 73°. The increase 
of the reflection component is numerically impor- 
tant (2.33 against 1.47 previously) but statistically 
unconvincing. The reflection is less predominant in 
the two other Seyfert classes. 

(b) Free inclination 

We then let the inclination free to avoid as much 
as possible the dependence between the different 
parameters. We find that Sy 1.5 galaxies would be 
seen closer to the normal to the disk with an in- 
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clination of 18° whereas the Sy 1 class would be 
viewed from an angle of 36° and the Sy 2 from 74° . 
It should be noted that in the PEXRAV implemen- 
tation of the reflection model, a stronger reflection 
hump systematically induces a large inclination. 

5.3. Investigating a link between the index slope F and 
the reflection parameter R 

We investigate whether the index slope T and the reflec- 
tion R are independently measured through the " draw" of 
this possible correlation using XSPEC (Figures 9 and 10). 



Confidence contour: 




Fig. 9. The 68%, 90% and 99% confidence contours of the 
index slope T and the Reflection parameter R for the Sy 
1 average spectrum. 



hard spectrum (e.g. r=1.6). We also note that if any class 
would have an index slope of less than 1.2, then reflection 
is not required to account for the data. 
The index slope is the signature of the primary X-ray 
source. Therefore, a link between the index slope and 
the reflection amount would indicate that the source 
of the cooling photons is related to the emission from 
the material at the origin of the Compton reflection. 
We will study using individual object spectra whether 
there is a correlation betwe en T and R in the Paper II 
i Deluit fc Courvoisierll2003|) . 



5.4. Results 

With a model independent method we find that Sy 2 
galaxies are similar to Sy 1.5 and differ from Sy 1 galax- 
ies. 

With a model dependent method, the best of all theo- 
retical models adopted for the two main classes is the 
PEXRAV model whereas the Sy 1.5 spectrum is well rep- 
resented by a broken power law. 

Considering the cutoff power law model, E~ r exp(-E/E c ), 
we find that Sy 1.5 and Sy 2 have similar behavior in the 
sense that their spectra do not require any cutoff whereas 
the Sy 1 spectrum presents a well constrained cutoff at 
70^7 keV. In model including reflection, the cutoff en- 
ergy in Sy 1 emission is found at 221_i5g keV and still 
not required in Sy 1.5 and Sy 2 emission. 
The main result is thus the presence of a high energy cut- 
off in Sy 1 that it is absent in other categories. The role 
and place of Sy 1.5 are not clearly identified. Indeed, Sy 
1.5 are close to Sy 1 concerning the relatively low reflec- 
tion influence, and similar to Sy 2 regarding the absence 
of a cutoff. 




5.5. Discussion and Perspectives 

Our results concerning the index slope of Sy 1 emission 
confi rm the result s obtai ned by IZdziarski et alJ l|l995h 
and iGondek et al.l 1^996) who found a~0.9. Our re- 
sults differ, however, in what regards S y 2 g alaxies. 
IZdziarski et all (|l995h and IZdziarski et"afl l|2000h found 
that Sy 2 are harder than Sy 1 which we do not find. 



Secondly, IZdziarski et alJ ( 1995h found a negligible reflec- 
tion in Sy 2 emission whereas in our study, we emphasize 
the possible importance of this component in the Sy 2 
average spectrum. 



Fig. 10. The 68%, 90% and 99% confidence contours of 
the index slope T and the Reflection parameter R for the 
Sy 2 average spectrum. 

This study is performed using an intermediate angle of 
cos 0=0.45 for all classes. For both Sy 1 and Sy 2 classes, 
if the average spectrum has a high index slope (e.g. T= 2) 
then it has a larger amount of reflection than if it has a 



5.5.1. High energy cutoff 

The presence or absence of a high energy cutoff is not in- 
fluenced by the inclination of the system with respect to 
the line of sight. The measured cutoff in the Sy 1 galaxies 
emission while none is found in Sy 2 galaxies clearly indi- 
cates intrinsic differences in their primary emission. 
However this lack could also be explained by a wide distri- 
bution of cutoff energies whereas the presence of a cutoff 
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in the Sy 1 average spectrum indicates a narrow range of 
comptonizing temperatures in the class. 
A cutoff in the Sy 1 emission provides, according to ther- 
mal Compton models, a reliable estimate of the electron 
temperature in the Comptonizing medium. The Sy 1 cut- 
off found at ksT=220 keV corresponds to a temperature 
of T~2.6T0 9 K for the electron gas in the case of unsatu- 
rated comptonization. 

If the absence of a cutoff in Sy 1.5 and Sy 2 galaxies is 
confirmed, it indicates either the importance of non ther- 
mal Compton scattering or a higher temperature of the 
Comptonizing medium. 

5.5.2. The role of reflection 

We show that Sy 2 galaxies may have a stronger reflection 
component than the other classes. Our fits also indicate 
that the inclination of Sy 2 could be larger than that of 
the other galaxies. Both results are not independent, only 
models free of this degeneracy will allow us to deconvolve 
the geometry of the reflection medium from inclination ef- 
fects. 

We also find that the average spectrum of Sy 2 galaxies, 
more affected by absorption than Sy 1, shows more re- 
flection than that of Seyfert 1 galaxies suggesting a link 
between the two processes. 

5.5.3. Implications of our results 

To investigate further these assertions, it is necessary to 
test our results in a different way, through the study of 
the distribution of the individual spectral parameters on 
a large sample of Seyfert galaxies. Therefore, we will study 
the distribution of the parameters T, E c and R for each 
object of a given class as far as can be measured to date 
and to compare them in order to test the results found in 
this paper. This study, presented in Paper II, will allow 
us to know whether the absence of a high energy cutoff 
in Sy 1.5 and Sy 2 is due to a wide distribution of values 
within the class and to study the correlation between the 
reflection and the power law slope. 

Seyfert 2 galaxies have different intrinsic emission and 
therefore different physical conditions in the inner regions. 
Sy 2 galaxies in which no cutoff has been found have prob- 
ably a higher temperature comptonizing medium than Sy 
1 galaxies. This reveals that the ratio heating/cooling of 
this medium is different in both classes. 
Our study has to be pursued within the Sy 2 galaxies 
where two subclasses known as Sy 2 with Polarized Broad 
Lines (PBL) and without PBL are found. Unified schemes 
predict that both types have a hidden Sy 1 nucleus. We 
have to search for common behavior in the intrinsic emis- 
sion of both Sy 2 subcla sses and to co mpare them with Sy 
1 emission features (see Deluit (2003)). 
The nature of the fueling mechanism occurring in the 
nucleus is unsolved and could also contribute to the 
Seyfert 2 particularity. Indeed, Sy 2 galaxies have been 



found to host starbur sts more often than type 1 AGNs 
llLevenson et allfeoOll) . We could expect that in Sy 2 ob- 
jects the accreted matter originates predominantly in star- 
bursts. 



6. Conclusion 

We obtained and studied the average hard X-ray spectra 
of Seyfert 1, 1.5 and 2 galaxies observed by the PDS on 
board BeppoSAX. We emphasized several differences in 
the central engine behavior of the two main Seyfert classes 
indicating that the Seyfert 1 and 2 galaxies have at least 
various physical conditions in their inner regions or are 
different objects. The differences cannot be due only to 
inclination/absorption effects as claimed by unified mod- 
els. 

The next steps in the Seyfert studies will be the search 
for the presence of a high energy cutoff in individual 
Seyfert galaxies to investigate the importance of the differ- 
ent physical processes in their emission (thermal and/or 
non thermal) and the investigation of the nature of the 
reflection medium (disk or material surrounding the nu- 
cleus), of its composition (density, homogeneity) and its 
localization. It is also important in a different set of con- 
siderations to investigate whether Seyfert 1 and Seyfert 
2 have intrinsic differences in their star formation rate, 
morphology and/or environment. 
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Appendix A: The final sample of the objects studied 

We present below the characteristics of the objects composing the final sample of our study. The Table A.l. includes 
the Seyfert 1 and Seyfert 1.5 galaxies subsamples and the Table A. 2. represents the Seyfert 2 subsample. 

Table A.l. General Characteristics and Spectral Properties of the objects composing the Sy 1 and Sy 1.5 samples 



O AT 

Source Name 
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T— 1 (j 
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^15-136 kev 
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Seyfert 1 £ 


;alaxies sample 






Mrk 335 


00 06 19.3 


+20 12 10.0 


0.025 


3.01cr 


1.38 


1.67 


Fairall 9 


01 23 45.6 


-58 48 14.0 


0.047 


11.8a 


5.42 


23.5 


NGC 985 


02 34 37.6 


—08 47 08.0 


0.043 


9.69a 


3.63 


13.1 


1H0419-577 


04 26 01.0 


-57 12 00.0 


0.104 


2.00a 


1.18 


25.7 


Mrk 110 2 


09 25 12.9 


+52 17 11.0 


0.035 


12.8a 


4.35 


10.4 
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10 34 38.8 


+39 38 34.0 


0.042 


2.14a 


0.71 


2.44 


NGC 3783 


11 39 01.7 


-37 44 20.0 


0.011 


48.5a 


14.2 


3.31 


NGC 4051 2 


12 03 09.5 


+44 31 50.9 


0.002 


7.25a 


2.06 


0.02 


NGC 4593 


12 39 39.4 


-05 20 38.0 


0.008 


23.0a 


9.63 


1.19 


MCG 6-30-15 


13 35 53.2 
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33.0a 


7.19 


0.84 
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21.9 


NGC 7469 
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28.0a 
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Seyfert 1.5 galaxies sample 


Mrk 1152 


01 13 50.2 
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3.32a 
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7.67 
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29.8a 
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Mrk 766/NGC4253 2 


12 18 26.2 
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14.0a 


3.63 


1.16 


NGC 5548 5 
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+25 08 12.0 


0.017 


52.6a 


10.6 


6.06 


Mrk 841 


15 04 01.2 
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2.71 


6.92 


NGC 7213 3 


22 09 16.2 
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0.006 


31.4a 


7.03 


0.49 



Note: 

1 Integrated Signal to Noise ratio in the PDS energy range 

2 2 available observations 

3 3 available observations 

4 4 available observations 

5 5 available observations 

6 Integrated Flux in 10 _11 ergs _1 cm~ 2 

7 Integrated Luminosity in 10 43 ergs _1 
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Table A. 2. General Characteristics and Spectral Properties of the objects composing the Sy 2 sample 
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01 


14 
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12 


25 
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012 


29.2i°J 


58.1a 


19.1 


5.12 


NGC 4939 


13 


04 


14 


5 


— 10 


20 


26 


9 
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18 


36 


58 


3 
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24 
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1.261 s ! 
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0.99 
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18 


38 


20 


3 
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25 
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7.96 


2.71 
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20 


23 


25 


6 
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31 


32 
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0.87 
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02 


02 


2 
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52 
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NGC 7314 
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0.30 
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23 
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21 


29 
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Note: 



Re ferences conce r ning the H ydrogen column dens ity: 

Bi lBassani et, alJ ljl999F) . Br: (Brandt, et alJ lll997fl. Ca : ICa"r>pi et alJ dl996h Co: IComastri et, all J199S 
Paper II o f this study, DC: iDella Ceca et al.1 <200ll), G: iGuainazzi et alJ jl998l). H: lHavashi et edJ r 
I2000J). M: lMaiolino et all (Il998l) S: ISmith fc Pond dl996l) . T: lTurner et al.l dl997h . U: IUeno et all ll99S 
1 19961) . W97: IWeaver et al.1 l)l997fl . X: IXue et, alJ <|l998l) 



LeiEhly^^d 



Weaver et al 



Hydrogen column density in 10 22 cm -2 units 

2 available observations 

3 available observations 
Integrated Signal to Noise Ratio 
Integrated Flux in 10 _11 ergs _1 cm -2 
Integrated Luminosity in 10 43 ergs _1 



Appendix B: The role of absorption in Seyfert 2 class studies 

According to actual models, the X-ray emission of many Seyfert 2 galaxies is characterized by a power law similar to 
that observed in Seyfert 1, with a cutoff at low energies (< 10 keV) due to absorption effects produced by gas column 
densities higher than 10 22 cm -2 . Se yfert 2 galaxies show a broad range of absorption, with column densities varying 
from < 10 22 cm~ 2 to > 10 25 cm -2 (Ris aliti et all fl999|) . These results yield to a further classification among Seyfert 
2 galaxies: Compton thin and Compton thick sources depending on whether their column density is respectively lower 
or larger than the inverse Thomson cross-section Njj= ov^, 1 =1.5T0 24 cm -2 . 

Therefore, for objects presenting a column density less than 1.5T0 24 cm -2 (Compton thin), high energy X-rays pen- 
etrate the absorbing material, making the nuclear source visible to the observer and allowing him to measure the 
Hydrogen column density. For objects with a column density higher than 1.5T0 24 cm~ 2 (Compton thick), the matter 
is optically thick to the Compton scattering and the nucleus is hardly visible. In the last few years, B eppoSAX data 
revea led a new highly absorbed class with only a lower limit on the column density Njj> 10 25 cm -2 ijRisaliti et all 

Il999l) . no detection of the intrinsic emission can be expected at energies between 15 and 15 keV. 

The p hotoelectric absorption cross section as a function of energy is well described by iMorrison fc McCammonl 
(1983) up to 10 keV. The component due to Compton scattering is described by the Thomson cross section equal to 
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6.65T0 25 cm 2 . Compton scattering plays an important rol e with the photoelectric cross se ction above 4 keV, and 
dominates above 10 keV. Starting with the analytic formulae iMorrison fc McCammonl 1 1983f) 



-24 — 2 

cm , 



^photoelectric = (C*0 + C\E + C*2-E 2 )-E 3 X 10 

where C =701.2, C x = 25.2 and C 2 =0. 
The photoelectric cross section at 10 keV is 

cr p/l (10keV) = 0.9532 • 10~ 24 cur 2 . 

At 10 keV the E -3 dependence is dominant, we thus extrapolate B.l. to higher energies using 

/ Energy > 



CTp/j (Energy) = cr p h(10keV) x 



V lOkeV 



(B.l) 



(B.2) 



(B.3) 



We present in the upper panel of the Figure B.l. the value of the total absorption cross section as a function of energy. 
The total absorption of the emission from the central engine is characterized as 



Flux* 



x g-orhJVff x e -<y P h{E)N H 



(B.4) 



We compute this at 15 keV, the lower energy used in our study, and show the result in the lower panel of the Fig. B.l. 
An object is considered as strongly absorbed when its flux has decreased by 50% due to absorption effects. We see in 
the lower panel of Fig. B.l. that this value corresponds to a column density around 70T0 22 cm~ 2 . 
Therefore, we consider that below a column density of 70T0 22 cm~ 2 the absorption is not an important effect at the 
energies covered in this study. Thus, we excluded objects with higher hydrogen column density. 
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Fig. B.l. Value and influence of the absorption at high energies. Fig. a gives the value of the total absorption cross 
section composed of the photoelectric and Thomson components. Fig. b gives the fraction of the rest flux at 15 keV 
as a function of the column density. 



Appendix C: Average Counts Spectra fitted by the different models 

We present below the counts spectrum of each class fitted by the different models applied in this study. 
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Fig. C.l. Sy 1 average spectrum fitted by a power law 
model and residuals 



Average spectrum of Sy 1 fitted by a broken power law 




Fig. C.2. Sy 1 average spectrum fitted by a broken power 
law model and residuals 



Average spectrum of Sy 1 fitted by a cutoff power law 
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Fig. C.4. Sy 1 average spectrum fitted by the PEXRAV 
model (cos 0=0.45) and residuals 



Average spectrum of Sy 1.5 fitted by a simple power law 
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Fig. C.5. Sy 1.5 average spectrum fitted by a power law 
model and residuals 



Average spectrum of Sy 1.5 fitted by a broken power law 
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Fig. C.6. Sy 1.5 average spectrum fitted by a broken 
power law model and residuals 



Fig. C.3. Sy 1 average spectrum fitted by a cutoff power 
law model and residuals 
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Fig. C.8. Sy 1.5 average spectrum fitted by the PEXRAV Fig. C.ll. Sy 2 average spectrum fitted by a cutoff power 
model (cos 0=0.45) and residuals l aw mo del and residuals 



Average spectrum of Sy 2 fitted by a simple power low Average spectrum of Sy 2 fitted by PEXRAV 




Fig. C.9. Sy 2 average spectrum fitted by a power law Fig. C.12. Sy 2 average spectrum fitted by the PEXRAV 
model and residuals model (cos 0=0.45) and residuals 



